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ABSTRACT 



Expressions for a generalized Modified 
Transitional Butter worrh-Chebyshev (HTBC) filter are 
derived. The characteristics of this filter as applied 
to digital filter design are investigated. It is shown 
that by adjusting location and order of the inserted 
zeros, the cut-off slope rate of the filter can be 
traded for naximuLn attenuation in the stop-band. 

The performance of this MTBC filter is compared to 
that of Butterworth, Chebyshev, transitional 
Butterworth-Chebyshev filters together with those 
suggested by other investigators C1]-C3]. It is shown 
that the step-band attenuation can be significantly 
increased without great sacrifice of cut-off slope 
rate. 

Step response of this MTBC filter is also obtained 
and compared with other filters. Various tabulations 
as well as graphs of this filter are given for design 
purposes. A computer program is developed for the 
design of this filter. 
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I. 



INTSODDCTION 



The digital filter is, as defined by Rabiner et al [10], 
'• a coaputational process or algorithm by which a digital 
signal or sequence of numbers (acting as input) is 
transformed into a second sequence of numbers termed the 
output digital signal". 

I 

The area of digital filtering can be divided into two 
major subdivisions as Finite Impulse Response (FIR) filters 
and Infinite Impulse Response (HR) filters. 

During the development of digital signal processing, the 
interest cf the investigators in HR and FIR filters varied. 
Before the introduction of the FFT algorithm by Cooley and 
Tukey (1965) HR filters were much more efficient than FIR 
filters. Stccham's work [13] on the FFT method of 
performing convolution indicated that implementation of 
high-order FIR filters could be made extremely 
computationally efficient ; thus, comparison between FIR and 
HR filters are no longer strongly biased toward the latter 
[5]. Because FIR filters require very high orders to 
produce a sharp attenuation shape, they are not often used 
for real-time filtering of waveforms. Recently, due to the 
increase in computing capabilities in digital signal 
processing and the availibity of long charge transfer device 
(CTD) tapped delay lines (TDL) , FIR filters are favored over 
HR filters. However, in applications like design of 
digital ccmb filters HR filters are the unique alternative. 

There are three basic design techniques of HR digital 
filters [5]. 



12 



First method is the direct design, which is, 
appropriately placing poles and zeros to approximate 
required frequency response. 

A second method is to use an optimization procedure to 
place the poles and zeros to match arbitrary frequency 
response specif ica tions. 

/ 

Finally the third technique makes use of highly advanced 
art of continuous filter design. This technique of 
designing digital filters from continuous filters by means 
of mathematical transformations is the most popular HR 
digital filter design technique. 

Standard Z-transform, Bilinear Z-transform and the 
matched Z-transform make possible direct transformation from 
S-domain to Z-domain, preserving essential characteristics 
of analog frequency response. 

Existence of frequency transformations reduces the 
problem to design a frequency normalized prototype low-pass 
filter. Then using appropriate frequency transformation, this 
prototype may be converted into desired band-pass, 
band-reject cr high-pass filter. Popular prototype filters 
are Burterworth, Chebyshev, elliptic and hybrid transitional 
filters. Frequency transformations for digital filters are 
discussed in various literatures ( [6] and [7] ). 

The problem of designing low-pass prototype filters, 
which possesses better stop-band attenuation and cut-off 
slope characteristics than existing prototypes has always 
attracted the researchers in the signal processing area. 

Budak and Aronhime suggested [1] modification of 
maximally flat rational functions by introducing a pair of 
finite transmission zeros such that the maximally flat 
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characteristic is maintained but the cut-off slope can be 
made steeper without great sacrifice of stop-band 
attenuation. 

Dutto Roy [2] investigated a more general case allowing 
insertion of multiple pairs of transmission zeros, either 
coincident or distinct. 

Introducing multiple pairs of jw-axis zeros in all pole 
Chebyshev transfer functions are investigated by Agarwal and 
Sedra [ 3 ]. 

The most attractive feature of these finite zero filters 
is that they offer the filter designer a great degree of 
freedom in choosing the location and order of the zeros to 
trade cut-off slope for stop-band attenuation. 

In this thesis, a modified Transitional 
Butterworth-Cbebys hev filter is developed, which is a more 
general case, introducing finite coincident or distinct 
multiple pairs of transmission zeros in transitional 
Butterworth-Chebyshev filrer. 

Trade-off's between the order of the filter, the order 
of transmission zeros, stop-band attenuation and cut-off 
slope are pointed out. Graphs helpful in the design of such 
filters are obtained. 

Performances of Butterworth, Chebyshev, Transitional 
Butterworth-Chebyshev filters and the designs suggested in 
references [2] and [3] are compared with those represented 
in this thesis for the orders of three through eleven. A 
computer program is developed to implement the filters 
mentioned above. 

In addition, the time-domain response of digital filters 
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is studied. There are many applications, such as digital 
ail filters, for which one is interested in the transient 
responses of filters that are specified in the frequency 
domain. Step responses of the filters that are discussed in 
this report are plotted and compared. 
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II. DE RIVATION OF UNIFIED TR ANSITION AL 
BDTTE RHOR TH -CHSBYSH EV FI LTERS 



A, INTRODUCTION 



The most popular technique for designing HR digital 
filters is to digitize an analog filter that satisfies the 
design specifications [5]. There are many techniques for 
designing analog low-pass prototype filters. Among the well 
known analog filter classes are the maximally flat 
(Butterwcrth) and equal ripple (Chebyshev) filters. 



Butterworth filters are simple, excellent in the 
pass-band and monotonic in both pass-band and stop-band. The 
Chebyshev filters are superior at and near cut-off frequency 
and at stop-band. 



The transitional Butterworth-Chebyshev (TBC) filters 
combine the desirable attributes of these two filters in a 
single approximation that is given by 




Z 



L 

it 




where 

k= Weighting factor 
n= Order of filter 

C .(w)= (n-k)— order chebyshev polynomial, 
when k=n, |F(jw) )2 is identical with the Butterworth 
function. When k=0, |F(jw ) |2 is identical with the Chebyshev 
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function. With a varying value of k, a TBC filter possesses 
seme characteristics of each. As k approaches to n, a TBC 
filter behaves more like a Butterworth filter, as k 
approaches to zero, it behaves more like a Chebyshev filter. 

In this chapter, modification of TBC filters by 
introducing finite coincident or distinct multiple pairs of 
transmission zeros will be discussed. It will be shown 
that, using a weighting factor k, and the location and order 
of inserted zeros as parameters, attenuation in the 
stop-band may be traded for sharpness of the cut-off 
characteristics. 

Expressions for the cut-off slope and minimum 
attenuation in the stop-band are derived in terms of order 
of the filter, weighting factor, location, and order of 
inserted zeros. 



B, MODIfICAIION OP THE TBC FILTE2S WITH COINCIDENT 
TBANSaiSSION ZEROS 



Introducing m identical pairs of transmission zeros at 
±jw<> to eg. (2-1) , we have 










In order to normalize |F(jw) |2, i.e. to force |F(jw) |2 to be 
egual to 1/2 at w= 1 , the constant K should be 
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Then eq. (2-1) becomes 






CWoVuD')‘"'+ Cu),‘.0^’" C„\(>o) 




where n>2m, because of the low-pass characteristics of the 
function, and the Chebyshev polynomial (w) is defined by 



Cos^ (n Cos"‘ui) > to4l 

Cosk^ (n &sV w) ^ , ui^l 



"• • Slope at cu t-off fre quenc y 



Substituting eq. (2-5) into eq . (2-4), we obtain 



Fcjui)! - 



Cu)o*- U3*’) 



2m 



<o3.'--ubV^"'+('^o"-0^"' CoS^[(n-k.)CoC'>£] 



Let 

h(u:>)= [(n-t) Cos'i^ 



(i-?) 

H 

H 



Using these values, eq . (2-7) becomes 



I ^ •{ L^) 
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we 



Taking the derivative of eg. (2-10) with respect 



get 






[-f '(to) K(ui) - 


W (10^ 


^ 1 1 [-f ^vo) (vOo^ - ^ 


V\(vjo)^ 



to 



w. 




At the cut-off frequency, w=1, we have 



f ’(w) = -^m (12.13) 

K ^U3) = i (2 -m) 

|F(3u,)| = i/lfT (i-lk) 



Substituting 
we obtain 

ko'^)! = 



these values into e 



no 

(u)o^-0 



4 - 



q. (2-11) and 

4 ( A - le.^ 



!L yo? 



simplif iying. 




For k=n , the result agrees with the cut-off slope of MB 
function, which is derived in [2]. For k=0, the result 
agrees with the cut-off slope of MC function as derived in 
[3]. 



3 • Stop -ban d ch aracte ri stic s 

In general, stop-band characteristics of finite zero 
filters will he of the form shown in FIG.1^. 

♦ Because of the inherent limitations in the plotting 
subroutine utilized to provide the graphs in this thesis, it 
was necessary to add supplementary axes to show proper 
scaling for seme of the graphs. 
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Figure 1 - GENERAL STOP-BAND CHARACTERISTICS OF THE FINITE 

ZERO FILTERS 
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The zeros introduced will cause a peak in the stop-'band, at 
a frequency Wp>w„ . The minumum attenuation in the stop-band 
may be defined as 

^ 0 i.0^ Tp 



where 




F 



In the stop-band 



eg. 



(2-4) becomes 




Let 

f(Of)= (uJo'-U3/)''"' 

W(u)f)= u/”" 





Taking derivatives of these values and substituting them 
into eg. (2-11) results in 



\ / 9 \2m 

I,-,. , 2 


2i\'1 / . t , t> 

lO ('0o-u3^J 




uJp JLn J 


-^1 






2CA-k)-2, 2ti 1 

-4 u3p J 


1^ 



(2.23^ 

Combining eg. (2-19) and eg. (2-23), we obtain 




n- ,2.m 





Thus 








(ir\-YV\-l2-0 



^ 



l(r\'2rv\) 

Ujo 
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And minumum stop-band attenuation will be given by 



C^M-rec = 10 n. m log (u)o) 



+»20 m ioj 



0 - SLrf\ 

n 



U3o^- 4 

LOo^ 







Plots of stop-band attenuation and cut-off slope of 
HTEC filters with two coincident transmission zeros, for 
orders 3 through 1 1 are given in figures 2 and 3. 



C. aODIPICAIION OF T3C FILT2BS WITH DISTINCT TfiANSMISSION 
ZZBCS 



Consider the n— order TBC function with m pairs of 
finite distinct zeros at±iwj , where i=1,2,...,m and w; >1. 



The aagnitude squared function of TBC function with 
distinct transmission zeros will be given by 



rw 






7 [ 



l*&\ 



i - 1 • - 1 



Let 



«vi 



3 M= T 

K.(u 3)= :f ( ui". 

I £ I 

Putting these values into eg. (2-27) results in 

“ i 



F(jw) 



i -+ 






H 
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Figure 2 - ST0P-3AND ATTENUATION OF MTBC FILTER WITH TWO 

COINCIDENT ZEROS 
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Figure 3 



COT-OFF SLOPS OF MTBC FILTER WITH TWO 

ZEROS 



COINCIDENT 



24 



I 


















1 • st o p -band and cu^of ^ characte r ist i cs 



occur, is given by 

= 0 






Id S'U)/ 



respect tc v, we obtain 

2 ,n~^ 

j-) ' 



__ a 



k^^io)[^U ^(>0) (i' 



solved for Wp may be found to be 



X COp 



m 




I 






• S.\ 



a ^ 0 



i-band peak 


will 




(±ii) 


eg. ( 2 - 28 ) 


with 


1 W (10) J 

K(to) i 


(2.30) 


eguation 


to be 




(.. 3 l) 


identical. 


Using 


attenuation 


peaks 



may be fcund to be 
■s=. o20 
where 

fft = |F(jWp;) 

and minimum stop-band atttenuation is given by 



min. = 



min. = rv\a"5^ 



5 -] 



(^. 2 ^ 



(2.2^ 
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Using eg. (2-5) and differentiating eg. (2-27) with respect 
to w, the cut-off slope becomes 




I 



+ Jl 





which agrees with eg. (2-17) when the W; *s are identical. 



Plots of cut-off slope and stop-band attenuation of 
MTBC filters with two distinct transmission zeros, for 
orders 3 through 11 are given in figures 3-17. 



D. SUHaARY 



The Transitional Butterworth-Chebyshev filter combine 
the best features of the Butterworth and Chebyshev filters. 
A Modified Transirional Butterworth-Chebyshev filter, 
obtained by introducing finite transmission zeros to a 
Transitional Butterworth-Chebyshev filter, possesses cut-off 
slope and stop-band attenuation, which are dependent on the 
modification parameters w© and m. The closer the w© is to 
unity, the steeper the cut-off slope. However, this 
improvement in the cut-off region results in degradation of 
stop-band attenuation. Thus, using w© and m as parameters, 
an advantageous trade between attenuation in stop-band and 
sharpness of the cut-off characteristics can be made. 
Graphs are given to serve as guides in trading cut-off slope 
for stop-band attenuation. 
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Figure 4 - CUT-OFF SLOPE OF MTBC FILTER WITH TWO DISTINCT 

ZEROS <w=1*06, n=5> 
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Figure 5 - COT-OFF SLOPE OF MTBC FILTER WITH TWO DISTINCT 

ZEROS =1#06^ n=6> 
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Figure 6 - CUT-OFF SLOPS OF J1T3C FILTER WITH TWO DISTINCT 

ZEROS < w =1*06, n=7) 
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Figure 7 - 



CUT-OFF SLOPE OF iilTBC FILTER WITH TWO DISTINCT 
ZEROS <w =1*06J, n=8> 
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Figure 8 - CUT-OFF SL0P2 OF MTBC FILTER WITH TWO DISTINCT 

ZEROS <w =1*06, n=9> 
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Figure 9 - CUT-OFF SLOPE OF MTBC FILTER WITH DISTINCT 
TRANSMISSION ZEROS <w =1*06 ,n=10> 
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Figure 10 - COT-OPF SLOPE OF HI 8C FILTER WITH DISIISCI 
TEABSalSSIOS ZEROS ‘ w ='»06 ,n-11> 
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Figure 11 - STOP-BAND ATTENUATION OF MTBC FILTERS WITH TWO 
DISTINCT TRANSMISSION ZEROS <w =1*06, n=5> 
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rigure 12 - STOP-BAND ATTENUATION OF MTBC FILTERS WITH TWO 
DISTINCT TRANSMISSION ZEROS <w =1-06, n=6 > 
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Jigure 13 - ST0P-5AND AITZNOATION 0? MT3C I1LI32S WITH T’liO 
DISTINCT THANS:iI3SI0N ZZHOS =1«06, n=7> 
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Figure 14 - STOP-BAND ATTENUATION OF MTBC FILTERS WITH TWO 
DISTINCT TRANSMISSION ZEROS <w =1*06, n=8> 
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rigure 15 - STOP-BAND ATTENOATION OF MTBC FILTERS WITH TWO 
DISTINCT TRANSMISSION ZEROS Cw =1#06, n=9> 
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Figure 16 - STOP-BAND ATTENUATION OF MTBC FILTERS WITH TWO 
DISTINCT TRANSMISSION ZEROS <w =1*06, n=10> 
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Figure 17 -- STOP-BAND ATTENUATION OF HTBC FILTERS WITH TWO 
DISTINCT THANSaiSSION ZEROS <w =1*06, n= 1 1 > 
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III. COM£ARISON OF MIBC FUNCTION W ITH B AND TBC 



A. INIBODUCIION 

In chapter 2 , derivation of MT3C filter_is given. It is 
shown that, stop-band attenuation and cut-off slope rate of 
this filter depend on modification parameters m, k, and Wo . 
So the filter designer have the flexibility to trade cut-off 
slope for stop-band attenuation by changing these 
parameters, without changing the order of the filter. 

In this chapter, performance of MTBC filter is compared 
with B, aB, C, ac, and TBC filters. 

Formulas for stop-band attenuations and cut-off slopes 
of all these filters are given in Table I. 

B. aiBC FUNCTION V.S. B FUNCTION 

Cut-off slope of HTBC function is given by 



FUNCTIONS 




(u)o*-0 



. 2 / ^ 
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From eg. (3-1) ratio of cut-off slopes of these two 
functions may be found as 



^TftC/S 



i^CvOo^-^) 



a 




v.s. w, for n=3 through 11 and m=1,2 are given in 
Figures 18 and 19. Stop-band attenuation of MTBC may be 
written as 






M76C 



0^0+ 4 ■+ c20 f(\ 




f\-2./Y\ 

a 



U)o^ 




Then the difference between stop-band attenuations is given 
by 



^M78c/a = ^ 






- 



— 4 (^n- rv\- +<^Orv\\o^ 



A-2m 

n 



W-4 

UJo^ 



(i-h) 

Plots of v.s. Wo for n=3 through 11 and m=l,2 

are given in FIG 20 and 21. 



Figures 19-21 indicate that, MTBC filter can be made 10 
times steeper than B filter still having 40 dB more 
attenuation at the stop-band. 



C. MTBC l.S. MB FUNCTION 



Using Table I, cut- off slope ratios 
filters may be found as 

' M7ac//M6 / , i" 

' 5m +a 4u)o--(; 



of MTBC and MB 




Plcts of w^ for n=3 through 11 and m = 1,2 are 

given in figures 22 and 23. 

The difference between the stop-band attenuations of MTBC 
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and MB is given by 




Equation (3-7) indicates that, stop-band attenuation 
difference between these filters doesn't depend on Wo . For 
n=10, m=1, k=l, with 48 dB more attenuation, cut-off slope 
ratio may be changed from 2.5 to 8.0 by changing w© from 
1.06 to 2.4. 

D. MTBC FUNCTION V.S. C FUNCTION 

The ratio of cut-off slope of MTBC function to cut-off 
slope of C function is given by 



Plot of Wo for n=3 through 11 and m=l ,2 are 
given in Figures 24 and 25. The difference between 
step-band attenuations of MTBC and C functions is given by 



given in Figures 26 and 27. The Chebyshev filter is known 
to provide much steeper cut-off slope than the corresponding 
B and TBC filters. Figure 25 shows that the cut-off slope 
of the MTBC filter can be made 1.4 times steeper than that 
of the Chebyshev filter, with n=5, m=2, and Wo=1.06. 

E. MTBC FUNCTION V.S. MC FUNCTION 





Plots of V.S w„ for n=3 through 11 and m=l,2 are 
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The ratio of cut-off slope of MTBC filter to cut-off 
slope of MC filter is given by 




Cu3o^--{) 

i ') 




Plots of v.s. Wo for n=3 through 11 and m=1,2 are 

given in Figures 28 and 29. 

The difference between stop-band attenuations of aiBC and MC 
functions is given by 

For a given Jc, is allways constant. 



F. MTBC V.S. TBC FOJ1CTIONS 



The ratio of the cut-off slopes of MTBC and TBC 
funcrions is given by 



'^T6c/T6C. = 1 + 



.m 



(u3o^-0 

Plots of v . s . w , for n=3 through 11 and 

given in Figures 30 and 31. 



(...) 



1=1 ,2 are 



The difference between the stop-band attenuations of 
these two functions is given by 



^ HT&c/t6C = — + <^0 



rt-2.rvi 



n 



Wo^ 



(3-1^ 



Plots of v.s. Wo for n=3 through 11 and m=1,2 

are given in Figures 32 and 33. 
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G. 



SOaMAfiY 



Using the location and the order of the inserted zeros 
and the weighting factor as parameters, the characteristic 
curve of the aTBC filter can be made steeper than that of 
the conventional all-pole filters without greatly 
sacrificing either stop-band attenuation or flatness in the 
pass band. A Modified Chebyshev filter provides slightly 
better performance than the MTBC filter, at the expense of 
substantial degradation of the pass-band flatness. Graphs 
are given to help in the comparison and in determining the 
numerical advantages gained by increased complexity. 
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Figure 18 - RATIO OF CUT-OFF SLOPES OF MTBC AND B FUNCTIONS 

V*S* W <M=1) 
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Figure 19 



RATIO OF COT-OFF SLOPES OF MTBC AND B FUNCTIONS 
7#S# W CM=2> 



49 



I 





{ 



i 







<2.t^eruio-\io^ ci'fferei^ce (^6) 




Figure 20 - THE DIFFERENCE BETWEEN THE STOP-BAND 
ATTENUATIONS OF HTBC AND B FUNCTIONS V*S* W <M=1> 
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Figure 21 THE DIFFERENCE BETWEEN THE STOP-BAND 
ATTENUATIONS OF MTBC AND B FUNCTIONS V*S« N <M=2> 
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Figure 22 - RATIO OF CUT-OFF SLOPES OF MTBC AND MB 
FUNCTIONS V#S^ W <M=1> 
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Figure 23 - RATIO OF CUT-OFF SLOPES OF MTBC AND MB 
FUNCTIONS V«S* W <M=2> 
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Figure 24 - RATIO OF THE CUT-OFF SLOPES OF MTBC AND C 

FUNCTIONS V*S* W <M=1> 
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Figure 25 



RATIO OF THE 
FUNCTIONS 



CUT-OFF SLOPES 
V*S* W <M=2) 



OF MTBC AND C 
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Figure 26 - THE DIFFE8ENCE BETWEEN THE STOP-BAND 
AITENOATIONS OF MTBC AND C FUNCTIONS V*S» W tM=1> 
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Figure 27 - THE DIFFESENCE BETWEEN THE STOP-BAND 
ATTENUATIONS OF MTBC AND C FUNCTIONS V»S» W <M=2> 
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Figure 28 - RATIO OF THE CUT-OFF SLOPES OF THE MTBC AND MC 

FUNCTIONS V*S« W < M= 1 > 
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Figure 29 - RATIO OF TH2 COT-OFF SLOPES OF THE MTBC AND MC 

FUNCTIONS V*S* W <M=2> 
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Figure 30 



RATIO OF THE COT-OFF SLOPES OF THE MTBC AND TBC 
FUNCTIONS V*S* W <M=1^ 
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Figure 31 - HATIO 0? THE COT-OFF SLOPES OF THE MTBC AND TBC 

FUNCTIONS V*S* W <a=2> 
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Figure 32 - THE DIFFERENCE BETWEEN THE STOP-BAND 
ATTENUATIONS OF MTBC AND TBC FUNCTIONS V«S* W <H=1> 
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IV. COMPUTER PROGRAM 



A. INTRODUCTION 

in the previous chapters, expressions for MTBC filter 
are derived and the performance of this filter is compared 
with various filters. 

In this chapter, A computer program is developed to 
implement MTEC filter. It is pointed out that all five 
filters, which are used to compare with the MTBC filter, are 
the special cases of the MTBC filter. Thus the program 
developed in this chapter for the implementation of the MTBC 
filter may also be used to implement any one of these five 
filters. 

A sample problem is worked out to illustrate the use of 
the plots presented in earlier chapters, the use of the 
computer pregram, and to point out the flexibility offered 
by the MTBC filter to the filter designer. 

B. COMPUTER PROGRAM 

In the most general form the transfer function of MTBC 
filter is given by eg. (2-27) , which is repeated here for 
convenience . 




where ; 



. 4Il . 

w = Location of i— inserted zero 

m = Order of inserted zeros 

n = Order of filter 

k = Weighting factor of Transitional Butterworth 
Chebyshev filters 

When k=n, and m=0, |F(jw ) |2 is identical to the 
Butterworth function. With k=n and m#0, |F(jw ) |2 is 
identical to the modified Butterworth function which is 
discussed by BUDAK and ROY [1] - [2]. When k=0 and m=0, 
|F(jw) 1 2 is identical to the Chebyshev function. With k=0 
and m#0 Modified Chebyshev function may be implemented which 
is discussed by AGARVAL and SEDRA [3]. 

When m=0 and k#0, |F(jw) |2 is identical to the 
Transitional Butter worth-Chebyshev filter. Finally allowing 
both k and m to vary. Modified Transitional 
Butter worrh-Chebyshev filters may be implemented. 

The program consists of two main parts. In the first 
part the analog filter, which is given by eg. (2-27), is 
implemented and the frequency response is plotted. 

In the second part, the analog transfer function F(jw) 
is first predistorted then transformed into z-domain by 
algebraic substitution method (using Bilinear 

z-transf ormation) to obtain H (z> . Digital transfer function 
[ H (z) ] is then factored into second order cascaded stages. 
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V . 3 . 



w, 20 



Finally, frecuency response curves ( I r ( jw) | 
log l?(j'^)l Wo and ?(jw) v.s. v© ) are drawn. 

Bilinear Transf oraa tion is preferred over rhe other 
available algebraic substitution nethods (i.e. Iipulse 
invariant, hatched z-transf ornation) in obtaining H (z) , 
nainly for the following reasons [4] ; 

(1) it has the property that realizable stable 
continuous systens are napped to realizable stable digital 
filters. 

(2) 'n'ideband sharp cut-off continuous filters can be 
napped to wideband sharp cut-off digital filters without the 
aliasing in the frecuency response. 

(3) After Bilinear 2r ansfornation, the relation between 
the analog and digital frequencies is given by 




3y choosing cut-off frequencies approaching to fj/2 where fi 
stands for sanpling frequency, extrenely sharp cut-off 
slopes nay be obtained. As an exanple plot of cut-off slo 
of bilinearly transforned Butterworth filter v.s. cut-o 
frequency is given in FIS. 34. 

C. 2BQUIRZD DATA CARDS 



The data cards required to use the program are given 
below . 

Card 1 : Values of n,a,!c in 312 format 
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Figure 34 - CUT-OFF SLOPE OF BILINEARLY TRANSFORMED 

BUTTERMOETH FILTER 
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Card 2 : Location of inserted zeros, w; where i=1,2,...,m 
in 8F10.5 format 

Card 3 ; Initial and final values of the frequency to be 
used for the frequency response plot in 2F10.5 format 

Card 4 : Number of solutions required in I3 format 

D. REQUIRED SUBROOTINES/FDNCTIONS 

In addition to the built-in subroutines, the following 
IBM source library subroutines are used. 

1 . POLRT 

2. PLOTP 

3. PSUB 

4. PMPY 

5. PADDM 

E. DESIGN EXAMPLE 



Suppose we want to design a digital filter with a cut 
off slope of 15 and minumum attenuation in the stop band of 
60 Db . 

Possible solutions for various types of filters are 
given in table II. 
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TABLS II* POSSIBLE SOLUTIONS OF DESIGN EXAMPLE 
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The following input data specifies a MTBC filter of 
order 7 with ^=2, Wo=1.36, k=1, sampling period of 1 sec., 
100 solution points, and frequency response plot from 0 to 3 
Hz. 



1 


M 5 


6 


7 1 8 1 9 1 10 


II 1 12 1 13 1 14 1 IS 1 16 j 17 t 18 1 19 1 20 


21 1 22123 124 1251 26 1271 28 |29| 30] 




7. 


i 


1 ! 1 


I 1 1 1 1 I 1 1 i 


t 1 1 t 1 t 1 t i 




1 ! 1 






1 1 1 1 1 \l\.\3t6 


1 i 1 ( 1 1 ) 1 t 




1 1 I 




o 


1 I I 1 1 1 i3i. lO 


1 1 1 1 1 1 1 i 1 


i 


o 

o 




111 


1 1 1 1 1 1 1 1 1 


L 1_[_. 1 1 1 1 1 1 




i 1 1 




1 1 L_ 


1 1 » 1 1 1 1 1 1 


1 1 1 t 1 1 1 t J 



Frequency responses of three of the possible solutions 
of the design example are given in Fig. 35. 
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Figure 35 - MAGNITUDE RESPONSES OF (a) MB FILTER N=8, M=1, 
(b) MC FILTER N=7, M=1, (c) MB FILTER N=7, M=2 
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V- T IM E DOMAIN RESPONSE OF DIGITAL FIL TERS 



A. INTRODDCTION 

In previous chapters we have investigated the problem of 
desining digital filters in the frequency domain, i.e. to 



designer should always consider the transient response 
characteristics of its filter. There are many applications, 
such as digital MTI filter, for which one is interested in 
the transient responses of filters that are specified in the 



characteristics of a filter will work against each other. 
Filters close to the ideal frequency characteristic can be 
designed. Filters whose time characteristic is close to the 
ideal can also be designed, but filters close to both cannot 



In this chapter time-domain response of digital filters 
will be discussed and it will be shown that the location of 
the transfer functions poles has a profound effect on the 
transient response of the filter. 

B. TRANSFER FUNCTIONS' POLES AND TRANSIENT RESPONSE 

Given a transfer function of a digital filter in 
factored form 



meet given frequency domain specif ica tions. A filter 



frequency-domain 



Time-domain and freque ncy-' domain 



[ 12 ] 
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where z, ’s and pj *s are the zeros and poles, respectively, 
of the filter. In general, the input signal's Z-transform 
is of the form 



)( (A ^ 

(i-q,) 




where g;'s are the poles of the input function. The 
response of the system may be obtained from the transfer 
function relationship 



(2,-qr) 

For simplicity, we assume that all of the poles of the Y (z) 
are distinct, and making partial fraction expansion of Y (z) , 
we obtain 



i(i) = H (4) •)((!%) = K. 



(i-p,)... 






4... + i!L±- 
i-pi 4-fn. 



The response of the system may be decomposed into [5] two 
parts called the input signal mode (ISM) and the system mode 
(SM) as 



V,:„ ii) = 



-+-••• + — 

i-qt 



/a) - -i- 

smU;- *r 

combining eg. (5-4) with eg. (5-5) and eg. (5-6) 



(s-c) 



i(i) = k, + y.sM (i) + Y sm ii) (s--?) 

Eg. (5-7) indicates that the response of any linear system 
to any input will contain modes generated by the input 
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signals poles and the system transfer function poles. The 
time domain response of the filter may then be found by 
taking the inverse Z-transform of eg. (5-7) 



Sguations (5-9) and (5-10) indicate that a fast responding 
system is one in which all of the system transfer function 
poles. Pi , are sufficiently smaller than unity in magnitude, 
in order that the system mode will decay to zero rapidly. 
On the other hand, a slowly responding system is one in 
which the system mode decays to zero very slowly (i. e. at 
least one of xhe p is close to unity in magnitude) . 

Those poles closest to the unit circle will be called 
the system's dominant poles because they tend to dominate 
the characteristic of the resultant transient response. 

The dominant poles' response property is, as stated by 
Cadzow [14] "The response time of a linear discrete system 
is directly dependent on the locations of the system 
transfer function's dominant poles. Depending on the 
particular response- time reguirement for a given 
application, we then have to correspondingly locate the 
dominant poles of the transfer function. A fast responding 
system necessitates dominant poles of magnitude much less 
than one ". 

The poles of the MTBC filter for the orders 3-11 and for 
various values of m and w^ are given in Table III. 





where 
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TABLE III* POLES OF !1TBC FILTER 
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TABLE III* CONTINOED 
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Although this partial fraction expansion method helps us 
to understand the importance of the poles of the systems 
transfer function in transient response analysis of digital 
filters, evaluation of residues of corresponding poles of 
partial fraction expansion is not a trivial problem. 

We believe that the so called ' transfer matrix ' 
method, which we are about to discuss, is more suitable for 
digital computer simulation. 



Given a digital filter weighting sequence h (n) 
input sequence x(n), the response of the digital 
be obtained by convolution summation ' 

rt. A. 

K-O Vtsro 



^ xCk.) =. ^ kCte.) x(n-k) 



and the 
filter may 




Eq. (5-11) may be written in matrix form as 
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where G is the systems transfer matrix. 



(f-ti) 

which is defined as 
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find 


the 


systems 


response 


using eq. 


(5-13) , systems 



transfer matrix must be available. To obtain the weighing 
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segueace (impulse response) of the system to form the system 
transfer matrix in terms of the transfer function’s 

coefficients, the difference equation of the system is 
solved for the impulse input, i.e.. 




m 

jCn-l) 

»=l 





where aj and bj are the numerator and denominator 
coefficients, respectively, of the transfer function of the 
digital filter. 

A computer program ( FORTRAN ) is developed to 
investigate the transient response of the digital filters, 
using equations (5-13) and (5-15). Program listing is given 
in Appendix B. 

Step responses of MTBC filter for various values of n, • 
m, and Wq are given in figures 36-39. Figure 36 and 37 
indicate that increasing rhe values of m and Wq also 
increases the overshoot and settling time, but doesnt have 
any significant effect on the rise time. The rise time tends 
to increase with increasing order of the filter. 
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c. suaaAEY 



Transienr responses of digital filters depend on the 
position of the poles of its transfer function. A fast 
responding filter has poles of magnitude much less than one. 
Hodification of all pole filters increase settling time, 
decrease peak overshoot, and doesnt significantly affect the 
rise time of the filter's step response. 
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Figure 36 - STEP HESPONSE OF ilTBC FILTER (N=7, W = 




1*46) 
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Figure 37 - STEP RESPONSE OF TNBC FILTER { N = 5, M= 1 , K=l) 
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Figure 38 - STEP RESPONSE OF MTBC FILTER ( M=1, W =1*36) 
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- STEP RESPONSES OF B, MB, 
a=1, W =1*36) 



Figure 39 



MTBC FILTERS ( N=5 



C /I ( T: 
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APPENDIX A 



C.OMPUTEH PROGRAM LISTING 



CIMENS ICN FNUMR ( 25 ) , RNUMI ( 25 I »RDNUMR (25 > ,RCMMI ( 25) . 
^Cl(25) ,C (25) ,Ri(25) ,F (25) 

CIMENSICN X(50),Z(50)»W(50) ,Y(50) 

CIyE^SIC^ FCCTR(50) t FCOT I ( 5C ) t COF( 50 ) 

CINcNSICN C^(50) ,CV(5C) »CZ(50) 

CIMENSICN RGQTDR(50) ,ROOTCI(50) 

CIMENSICN CYK50) 

CIMENSICN FM 100) tPV ( ICC) tFC (ICO) ,RC1( 100) 

CIMENSICN FN(25 ) ,RC ( 25) 

CIMENSICN >T (25) ,YT ( 25) ,AX (25) »AY(25 ) 

CIMENSICN CUMMY{ 100) t DUMMX ( 100 ) , I TB ( 23 ) »RTB(30) 
CIMENSICN CLMKAY(ICO) 

CCMPLE) F1,R,C,C1 

REaC the CFCEF CF filter, GRCER CF inserted ZERCS AMD THE 
WEIGHTING FACTCR 

F£AC(5,5C1) N,M,K 
5C1 FORMAT (212) 

FEAC LCCATICN CF ZEROS 

F£AC(5,5CC) (rt(I),I=l,M) 

5CC FORMAT (EFIC .5 ) 

RE/0 INITIAL ANC FINAL VALUES CF THE FRECUENCY 
(10 BE USED IN FREQUENCY RcSFCNSE CALCULAT ICNS ) 

READ(5,5C2) WBEGIN,WLAST 
5C2 FORMAT (2F1C. 5) 

RcAC THE NLMEEP OF SOLUTION POINTS 

R£AD(5,5CE) NPCINT 
5CE FCRMAT(I2) 

REAC SAMPLING PERIOD CF THE DIGITAL FILTER 

READ (5,70 T 
7C FORMAT (F 1C. 5) 

XPCI NT=NFCINT 

WCELTA= (WLAST-WBEGIN ) /XPOINT 



CALCLLATICN CF NUMERATCF FCLYNCMIAL 



IDIMY = 1 
V ( 1) = 1 .C 

IF (M .EC .C ) GO TO 401 

CO 2 1 = 1, M 

X(1)=W(I )X=»4.0 

X ( 3 ) =-2 .C^W ( I )«*2.C 

X(2) =C.C 

X(5)=1.C 

X(4)=0.C 

CALL FNFV (2,IDIMZ,Y,ICIMY,X,5) 

CO 3 J= 1 , IC IMZ 
V ( J) =Z ( J ) 

2 CGNT I NLE 
ICIMY=ICIMZ 
2 CCNTINLE 

CCNVERSICN OF INDEPENDENT VARIABLE FROM 



CO 26 1 = 2 ,ICIMY,4 
Y (I )=-l.C^Y (I ) 



W TO S 
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2t CONTINUE 

FIND THE FCCTS OF NUf'ER^TCR S^U^REO FUNCTION 
NNUM=I CIKV-1 

CALL PCLRT ( V ,C0 F , MNL , RO CTR , RCOTI , I ER ) 



SELECT THE RIGHT HALF PLANE RCOTS OF NUMERATOR 
FUNCTION 



L = 1 

CO 20 I = 1,MNU.M 

IFCRCCIF (I I.GT.O.O) GO TO 20 
RNU.MR ( L ) =FCGTR( I ) 

RNUMI ( L > = RCCTI ( I > 

L = L+1 

2C CONTINUE 
VN=MNUM/2 
LK = L 
LZ=L-1 

IF(LZ .EC ) GO TO 2 2 
LL = L .0 

CO 35 I=1,MNUM 

LX=L+1 

LT=I+1 

IF (RNUN I (LL) .EQ.ROQTK I ) .OR .RNUM I ( LL > . EG . ( - 
’(I))) GC 7C 35 
IF (LL .EC. 1) GO TO 26 
LJ=LL-2 

IF (RNUM (LJ) .5G.RCCTI ( ! ) .QR.RNUMKLJ ).EG.(- 
=♦(!))) GC TC 25 
IF(LL .EC .2 ) GO TO 26 
LI=LL-^ 

IF (RNUM I (LI ) .5Q.RCCTI (I ) .0 F .RNUM I ( LI ).EG.(- 
»(I))) GC TC 25 
26 RNUMR ( L ) =FCCTR(I ) 

RNUMI (L ) = RCCTI(I ) 

RNUMR(LX) =FCCTR(LY ) 

RNUMI ( LX ) =RCCTI ( LY ) 

IF(LX .EC .MN ) GO TO 22 
LL=LL+2 
2: CONTINUE 
32 CC 22 1 = 1, MN 

R(I)=CMFLX(FNUMR(I ) ,PNUM (I) ) 

22 CONTINUE 

CALL MAKFCL (MN,R,C) 

MM = MN + 1 

C(MN1 ) =CMFLX( 1.0, 0.0 ) 

CG 31 1 = 1, MM 
RC(I )=REAL(C(I) ) 

21 CONTINUE 
GC TC 5 
^C1 RC(1)=1.C 
MM=1 
MN=C 

RNUMR ( 1 ) =C .C 
RNUMI ( 1 ) = 0 .0 

CAlCULATICN OF CENOMINATCR POLYNOMIAL 



e 



1 1 
7 



NDNUM = 2=»N 

IF (M .EC .C .ANC.K.EC .N ) GO TC 402 
CCNST=1 .C 

IF(M.EC.C) GO TO 7 

CC 11 1 = 1, M , ^ ^ 

CMULT = (V ( I )=»*2 .0-1 .0 J=»*2.0 

CONST=CMULT*CONST 

CONTINUE 

IF(K.EC.C) GC TO 52 



KK = 2=»K 
KKK = KK ^ 1 



SCUARED 



.O-ROOTI 

l.C=*'ROGTI 

.0=*RQCTI 
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Cf«‘(KKK )=CCNST 
CG 12 I=ltKK 
C*V(I )=C.C 
12 CONTINIE 
GG TC 52 
52 CM( 1 ) = 1.C 
KK = 0 .C 
KKK=1 ,C 
52 M<=N-K 
NKK1 = NK + 1 

OBTAIN THE CHEBYSHEV FCLYNCMIAL tSGUARE IT AND HULTIPLY 
THE RESULT BY BUTTERHORTh SQUARED FUNCTION 

CALL CFES V (CYl ,NK) 

CALL FNPY (CY,NKK ,CY1,NKK1 ,CY1 ,NKK1) 

CALL FNFY ( C2 , 1 DI MC Z » CY ,NKK ,CM , KKK ) 

CETAIN CENCHINATOR SQUARED FUNCTION 

CALL PACCN (Z,IDIHZtY,IDIMY,1.0,CZ»IDINCZ) 

CCNVERSICN CF INDEPENDENT VARIABLE FROM W TC S 

^ CG 17 I=2,ICIMZ,4 
Z(I )=-l .C^Z( I ) 

17 CGNTINLE 
£ NDNUN=ICIN 2-1 

FIND THE FCCTS OF DENOMINATOR SQUARED FUNCTION 

CALL PCLRT ( Z,COF , NONUM , ROOTOR ,ROOTDI » lER ) 

GC TC ^C2 

4C2 CALL RCCT (NCNUM, RCCTCP tPCOTCI ) 



SELECT THE RIGHT 
FLNCTION 



HALF PLANE ROOTS OF DENOMINATOR SQUARED 



ACE J=i 

CO 21 I=1,NCNUM 

IF (RCCTCP ( I ) .GT .0 .0 ) GO TO 21 

RDNUMR ( J ) =FCCTOR ( I ) 

RONUMI(J)=RCCTDI(I) 

J = Jfl 

21 CONTINUE 
MC=NCNLM/2 
DC 23 1 = 1, NC 

RKI ) = CMFLMRDNUMR(I ) ,RDNUM (I )) 

22 CONTINUE 

CALL HAKFCL (MC,R1,C1) 

HC1 = MC-»1 

CKMCl >=CHFLX(KC,C.C > 

AT FORMAT (2> ,25F11. 3// ) 

CO 32 1 = 1, MCI 
PCK I )=P£AL(C1(I ) ) 

32 CONTINUE 

NCRMALI ZAT ICN CF NUMERATOR POLYNOMIAL 

FACT0R=RC1( D/RC ( 1 ) 

CC ICC 1=1, MM 
RC (I ) =FC (I )>>FACTOR 
IOC CONTINUE 

R I T E ( £ , 5 7 C ) 

OUTPUT SECTION 

57C FORMAT (’ 1 • ////) 

5£C FORMAT (2X,‘ LOW-PASS PROTOTYPE (CONT I NUOUS ) FILTER'//) 
HRITE(£ ,571) 

571 F0RMAT(2X,' CRDER CF FILTER '//) 
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6CC 



573 

574 
SCI 
576 



CF INSEkTSO zeros '//) 



>RIT£(6,gCC) N 
FSRM/T (2X , 12// ) 

URIT5(6 , 57 2 ) 

572 FORMAT (2> , ' ORDER 
V^PITE (6 ,£CC ) M 
V»RITE (6, 573 ) 

FORMAT (2>,' WEIGHTING FACTOR OF TBC FILTER •//) 
WRIT£(6,SCC) K 
WRITE( 6,574) 

FORMAT (2>,‘ LOCATIONS OF ZEROS '//) 

WRIT5(6,£C1) (w(I),I=l,N) 

FORMAT (SFIC. 5//) 

WRITE ( 6 , 57 6 ) 

FGRMAT(2X,’ COEFFICIENTS CF NUMERATOR PCLYNCMIAL 
=»( CESCENCING ORDER ) '//) 

WRITE(6,47) (RC( I ) , I = 1,MM ) 

WRITE( 6 , 57 5 ) 

575 FQRMAT(2X,' COEFFICIENTS OF DENOMINATOR POLYNOMIAL 
» ( CESC6NC ING CROER ) ' // ) 

WRITE(6,47) (RCKI ) ,I = 1,MC1) 

WRITE (6 , 51 1 ) 

51] FORMAT (2X , 'ROOTS CF NUMERATOR'//) 

«RIT£(6,512) 

512 FORMAT (6X ,' REAL P ART ' , AX , • I NAG IN ARY PART '//) 

CO 52C 1 = 1, NN 

WRITE (6 ,512 ) RNUMR( I ) ,RNUVI ( I ) 

513 FORMAT (5X ,F10.5,4X ,F1C.5/ ) 

52C CONTINLE 

WRIT£(6,514) 

514 FORMAT (2X , 'ROOTS CF CENCNI NATGR' // ) 

WRITE (6, 5 12) 

CO 521 I=1,NC 

WRIT£(6,513) RONUNR(I) ,RDNLNI(I) 

CONTINLE 
WRITE(6,576) 

FQRMAT(2X,' INITIAL VALLE CF FRECUENCY '//) 
WRITE(6,51C) WBEGIN 
51C FORMAT (2X .F1C.5//) 

WRITE ( 6 , 57 7 ) 

=CRMAT(2X,' FINAL VALUE OF FRECUENCY '//) 
WRITE(6,51C) WLAST 



521 

576 
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PLOTTING SECTICN 



WRITE ( 6 ,550 
55 C FGRNAT('l') 



CALL PPLCT TO 
F ILTER 



FIND FRECLENCY RESPONSE VALLES CF THE 



CALL PFLCT (RC ,RC 1, W5EGIN,'WCELTA,NP0INT ,PX, PY ,ND1 , NNl ) 
CALL FLCTF ( ?X , PY , NP C INT , C ) 

PRECISTGRTICN FOR BILINEAR TRANSFORMATION 
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SCALE=1AN (T/2. ) 

CO 777 LK=1,NN1 
£XF=LK-1 

RN(LX)=RC(LK)/(SCAL£-*EXP) 
CONT INLE 



CO 776 LK=1,ND1 
EXF = Li<-l 

RD(LK) =RC1(LK)/(SCAL£=*^EXF) 

776 CONTINLE 
NT = MN 
M K = M0 

CALL ZCNN (RN, RO, N K, NT, WBEGIN, WLAST, WCELTA,NFCI NT, 7) 
CALL XFCFN(RC,RC1,1,XT,YT,NN,NC,100C«C,12»0 
CALL FFLCT(XT,YT,5C0.C,10.C,100,AX,AY,NC1,NM) 

CALL FLCTF (AX, AY , ICO ,0 
STOP 



88 



OOO OOOOc^OOOOOO 



ENC 



c 

c 

c 

c 

c 

c 



EL3RQI.TINE CH6SV (CY,NC) 

O-cEVShEV 

NC : CRCEF Cf Ch£EYSI-EV FCLYNCMIAL 

CY : CALCILATEC COEFFICIENTS (CESCcNOING CRCER) 

CI/ENSICN CX(50) ,CY(5C),YY(50),ZZ(50),Z(50) 
NN=NC + 1 
CX(1 ) = 1 ,C 

IF (NC .EC .C) GO TC 2 
CY (1)=C.C 
CY(2 > = 1.C 

IF(NC.EC.l) GG TC 19 
YY ( 1 ) = C .C 
YY(2J=2.C 
CO 5 1=2, NC 
II=I+I 
11 = 1-1 

CALL PNFY ( 2 , IDIMZ , YY ,2tCY, I ) 

ICIMZZ=ICINZ 
CC 6 J=1,ICINZ 
ZZ(J ) =Z( J) 

E CG.NTINLE 

CALL FSL8(Z,IDI^<Z,ZZ,ICINZ2,CX,I1) 

LFCATE FCLYNCMALS 



IE 



CQ 7 J=1 , I 
CX ( J ) = CY ( J I 
7 CONTINLE 

CO 8 J=1,ICINZ 
CY ( J ) = Z ( J I 
£ CCNTINLE 
■ CONTINLE 
r ETURN 
CY ( 1 ) =CX (1 ) 

GC TC IE 
ENC 



SL'SRCUTINE ZCMN ( R C , R C 1 , N , V , F 1 , F2 , WOELT A , NP , T ) 
SLcPCUTINE TC FIND ANC FLCT F R EQUENCY -A E S PCNS = CF DIGITAL 

filters 

RCl : COEFFICIENTS OF CENCNINATCR CF CONTINUOUS FUNCTION 
RC COEFFICIENTS CF NUMERATOR CF CCNTINUCLS FUNCTION 
N : ORDER CF CENCMINATCR 
N : CRCER CF MNERATGR 
FI : INITIAL FRECUENCY 
F2 : final FRECLENCY 
7 ; SAMPLING PERIOD 

DINENSICN ECC(25) 

DIMENSION RC<25»,RC1(25»,ZX(25),ZX1(25),XJ(25),YJ(25), 
*RCGTPC(25) ,RCGTRN (25 ) ,RCCT IC(25 ) , ROOT IN (25 ) ,CY(25 ) , 

^ XPAY( 2 f ) , >FAYOA( 25 ) , 

= CX(25I,YFAY(3> ,YPAYDA (2 ) 

Clt*Et^SlCS ACC(iO) ,Ail(10),A22(10),Bll(10),c22(10) 

DIM ENSIGN XAXIS(15CO),YAXISi(15O0),YAXIS2(15CO), 
=*'FF(15CC ) ,CEL(150C) 

DIMENSION >CCF (3 ) , X2CCF(3 ) ,CGF (3 ) , RGCTR (21 , RCCTI ( 2 ) 

1 ,RCOTl (2 ) ,RCCT2(2 ) 
lIMENSICN CLNKAY(ICO) 

DOUBLE FRECISIGN 
CGMMCN/RAW/VsdOO ) ,Y( ICC) ,ICALL 
REAL M AGF ,.*>AGNF, NLMi ,NLM2 
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c 



c 

c 

c 

c 



CCMPL£> C>,CY,XJ,VJ 
CCf^PLEX MNtCENOM 
MJ = M+ 1 
NJ = .N + 1 
NV 

CALL E 2>FPy (RC »M ,Nyv , ZX , 1 ) 
CALL E2XFFN(RCl,N,Nf^iV,ZXl,C) 



9C 



SFACT=ZX(^J)/ZX1(^J) 

CC 9C 1 = 1, NJ 

ZXK I ) =ZX1 (D^SFACT 

CCNTINLE 



VPITE ( 6,60 



CCEFFICIENTS CF CENOMINATCR PCLYNOf-IAL 



VRITE (A ,45 nZXK n , I = 1,NJ ) 
VRITEIA,AC) 

CALL ‘ ‘ 



CALL 
CALL 
CALL 
CALL 
F = F1 
CO 1C 
X = F=*T 



PCLR7 (ZX.COF, N, FCOTRN ,RCOTIN, I £R ) 
FCL = T(ZX1,C0F,N , RCQT R C , RGCT 10 , i ER 
CCE VEiN (iN ,MX) 

FACTCR (ROOTRN,RCCTI.\ , N , MX , AOO , A 1 1 
FACTCR (RQQTRC ,RCCTIC, N.'^X ,322, Eli 



1 = 1 , Ni 



) 

, A22 » 
, ECO) 



MAGh = l ,C 
FF3 = C ,C 
CQ A J = 1 ,>X 
AC = ACC ( J ) 
A1=A 11 ( J ) 



A 2 =A 22 ( J ) 

El = 311 (J ) 

E2=E22 ( J ) 

NLM1 = ( A2 + AC)=«=CCS( X)+A1 
NUM2= ( AC-A2 J-SIN ( X ) 

NLV=CVFLX (^LM1,NUM2) 

CEN0M1=(E 2+1.0 )*CCS(X)+E1 
CENGN2 = ( 1 ,-E2 )-S IN ( X ) 

C5NG.'-=CMFLX(DENGV1,CENCM2) 

MAGH =A ES ( ^AGh’^CA 3 S (NLM/OENGM) ) 

FH1= (A7AN (^LV2/NUM1 ))^57.2957 7S5 1 
?H2= (ATAN (CENCM2/CENCM) ) *57 .29 57795 1 
Fh2=FF '♦ ( FF1-PH2 ) 

A CCNTINLE 

MAG.Nh = 2C .+ALCG10 ( VAGF ) 

Fh ( I ) = FF2 
FK = Ft- I I ) 

>AXI S ( I ) = F 
YAXISI ( I )=NAGH 
YAXI S2 ( I ) =NAC-NH 
(I ) = X A X I E ( I ) 

Y(I)=YAXIS1(I) 

F=F+VC£LTA 
1C CCNTINLE 
■f.PITE (A ,AC ) 

AC FORMAT ( ' 1 ' ) 

45 FORMAT (2X , 10F12. 2/// ) 

WRITE (A ,7C ) 

7C F0RMAT(2X,' BILINEARLY TRANSFORMED LCW-PASS 
*(FRCTCTYFE) DIGITAL FILTER'//) 

WRITE ( A ,A 1 ) 

A1 FCRMAT(2X,' COEFFICIENTS OF NLMERATOR FCLYNCMIAL 
WPI7E(A,45 MZXd ) ,I=i,NJ) 

WRITE ( A , A2 ) 

A2 FORMAT (2X ,' COEFF ICIENTS OF OENGMINATCR FGLYNCMIAL 
WRITE ( A ,45 ) (ZXl ( I ) , I = 1 ,NJ ) 

WRITE ( A , AE ) 

A3 FORMAT (2X,‘ NUMBER CF CASCADED STAGES’//) 
WRITE(A,A4) MX 
A4 FORMAT ( IX , 13 ) 



' // ) 



'//) 
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66 
S6 
c c 

12 

51 

16 
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lA 



VSITE(6 ,6 1) 

(2X, ^COEFFICIENTS OF CASCADED STAGES'//) 

»^22(I ),eil( I) ,822 ( 1 ) 

FUS'^fiT (iX,5rl2.5//) 

CENTIME 
’.^RIT£( 6 ,65 ) 

FORMAT (2X , 'FCGTS CF NUMERATOR'//) 

CC 66 I'=1,N 

ViR IT E ( 6 , ^5 ) ROOTRN ( I ) , ROOT IN ( I ) 

CCNTINL E 
I'RITE ( 6 ,67) 

FCRMAT(2X,' ROOTS OF C5NOMINATCR '//) 

CC 66 1 = 1 ,N 

i^SITEC 6 ,^5 ) ROOTRDd ) , ROOT ID (I ) 

CONT INL E 
WRITE (6 ,S6 ) 

FORMAT(2>,' sampling FRECUENCY '//) 

WRITE(6,SS )T 
FORMAT (2X ,F1C.5//) 

N=N? 

W R I T E ( 6 , 1 2 ) 

FCR*^ AT Cl') 

CALL FLCTF OAXIS ,YAXIS1,N,C) 

W R I I E ( 6 , 5 1 ) 

FORMAT ( //ACX , ' A = S. GAIN V.S. FRECUENCY') 

W R I T E ( 6 , H ) 

FORMAT (' i ' ) 

CALL FLCTF ( XAXIS,VAXIS2,iN,C) 

WRIT£(6,52) 

FORMAT ( //4CX , ' GAIN (C3 ) V.S. FREQUENCY ') 
WRITE ( 6 , lA ) 

FORMAT ( ' 1 ' ) 

CALL FLCT F ( XAXIS , Fh ,N ,C) 

WRITE( c ,52 ) 

FORMAT (//ACX , ' PHASE V.S. FRECUENCY ') 

RETURN 

ENC 



SLERCUTINE CCEVEN(N,M) 

C SLnROU'TINE TC FIND THE NUMBER CF CASCAC6C STAGES 
C N: CcGREE CF NUMERATCR FCLYNOMIAL 
C M: degree CF CENCMINATCF FCLYNCMIAU 
C 

NN = (N/2 ) ^ 2 
IF (N .NE .NN ) GO TC 1 
M = N/2 
GC TC 2 

1 v=(n+1)/2 

2 RETURN 
ENC 



SUB ROUTINE £ ZXFRM ( X , M ,NMM , ZX , I FACT) 

DIMENSION >(25),Y(25),ZX(25) ,XCUM(2) 

’ , YCUM( 2 ) , 2 1 (25 ) , Z2 ( 25 ) ,X1 ( 25 ) , X2 ( 25) 

=* ,XFACT ( 25 ) , 2X1 (2 5 ) 

C SUBROUTINE TC FIND THE BILINEAR Z-TRANSFGRM OF A GIVEN 
C FClYNCMIAL 
C X: FCLYNCMAL 
C : CRCER CF FCLYNOMIAL 

C NMM=N-'^ { THE DIFFERENCE BETWEEN THE ORDERS OF NUMERATOR 
C AND DENCMINATCR 
C 2X: RESULTANT FCLYNOMIAL 
C IFaCT= 1 FCR NUMERATOR 
C IfACT= C FCR DENCMINATCR 
C 

^y = ^^ + ^^,v+l 
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N2=H+1 
DO 1 

ZX(I )=C.C 

1 CQNTIME 
XOt.v ( 1 ) =-l .C 
XDUy (2 )-=l.C 
V0L.V( 1 ) = 1 .C 
Y0LM(2 ) = 1.C 
IZX = 2 

CO 2 11 = 1, ^2 

CALL PCLE>F(XDUM,2,yN,Xl,I 1) 

CALL PCLEXF (VCUy ,2 ,X2 , 12 ) 

CALL F^PY( 21, IZ, XI, I 1,X2,I2) 

FACT = X (II ) 

CALL PADCM Z2 , 1 Z2 , ZX , I ZX, FACT, Z1 , 1 Z) 
CC 4 K=1,I22 
ZX(K ) = 22 (K ) 

4 C3NTINLE 
NM=Mf' + l 
N = «-l 
IZX=IZ2 

2 CCNTINLE 

IF(IFACT.EC.C) GC TC 5 
CALL PCLEXPCYDUM, 2,NyM,XFACT,iNXFACT) 
CALL P^FY<ZX1,MZ,XFACT,NXFACT,ZX,IZ2 ) 
DC 6 I = 1 2 

ZX( I ) = ZX1 ( I) 
i CCNTINLE 
f V=M2-1 
PETURN 
END 



SLeRCLTINE FACTOR ( R R N , P IN , M , VX , A22 , A 1 1 , AGO ) 

CIMENSIGN RRN(25) ,RIN(25) ,XR(25) ,XI (25) ,YR(25),YI (25) , 
C(25),C1(25 ),A00(25 ),A11(25),A22(25) 

CCN.PLE> C,C1 

INITIALIZE CCLNTERS 



11 

16 



16 



IM = C 
IY = l 
K = 0 
M=0 
IX = 1 
NXR = C 
NXI = C 

IF(M .NE .C ) GC TO 18 

CC 17 1 = 1, NX 

AOO( I )= UC 

A11(I)=C.C 

A22(I )=C.C 

CCNT INL E 

GC TC 2C 

IF(M .NE . 1 ) GC TO 16 

AOO(l»=l,C 

AIK 1) =-l.C*RRN( 1 ) 

A22( 1)=C.C 
K = 2 

NOIF = N»> 

NXK= 1 
GC TC 1C 
CO 1 I = 1 , ICC 
IX1 = IX-H 

IF(IX.GT.N) GO TC 7 

IF(RRN(IX).EC.RRN(IX1).AND.RIN(IX).EG. 
^ (-1.0=*FIN(IX1) )) GC TC 2 
IN=I N+1 

XR ( IN I =-l .C’^RRNl IX) 

IX=I X+1 
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NXR=NXP+ 1 
GC TC 1 

2 m=iv+i 

YR( IV )=PRN (IX > 

YI (IY) = FU (IX) 

YRdYl )=RRN(IX1) 

YI (lYl ) = RI^( IXl) 

IX=IXlHl 
M=NI + 2 
NXI=NXI+1 
IY=I Yl+1 
1 CONTINUE 
7 NT=NXP4NXI 

IF(NT.EC.NX) GO TC 3 
IF(MX.GT,NT ) GO TO S 
ncif=nt-n> 

IC = 1 

CO 4 I=1,NCIF 
IC1=IC->1 
ACG( I ) = 1,C 

^ 11 ( I )=XR ( IC )+XR ( ICl ) 

A22( I ) =XR ( IC )*XR ( ICl ) 

IC=IC+1 
^ CONTINUE 
NOIFl=NCIF+l 
NCGF=NCIF*2 
NUAST=NX-NXI 
CC 5 I=NCIF1 ,NUAST 
NC0F=NCCF4 1 
ACO( NC IFl ) = 1 .0 
AlKNCIFl )=XR(NCCF) 

A22(NCIF1 ) = C .0 

5 CONTINUE 
NeEGIN=NLAST+l 
NCOMF= I 

CO 6 I=NEEGIN,MX 
NC0MP1=NCCNF41 

C ( 1 ) =C^FLX (YR (NCONF ) ,YI (NCCMF ) ) 

D( 2) =C ^FUX (YR(NC0NF1 ) , Y i ( N CCM P 1 ) ) 
CAUL NAKPCL(2,D,D1) 

AOC(NEEG IN ) = 1.0 

AIK NBEGIN ) = FEAL (Cl (2) ) 

A22(NEEG IN ) = REAL (Cl ( 1 ) ) 

NCOMF = NCCN Fl + 1 

6 CONTINUE 
GC TC 2C 

f NCIF=NX-N7 

IF(NXR.EC.C) GO TO 15 
CC 9 1 = 1, NXR 
A00( I ) = 1 .C 
AIK I ) =XR ( I ) 

A22( I )=C ,C 
K = K+ 1 



c CONTINUE 

IF(NXI,EC,C) GO TC 1C 
15 NXK=NX? 

CO 11 1=1, NXI 
NXK = NXI< + 1 
IK=I+1 

C ( 1 ) =CNFUX (YR( I ) , YI ( I ) ) 
0(2)=CNPUX (YP(IK) ,YI (IK) ) 
CALL MAKFCL (2,0, Cl) 
ACC(NXK ) = 1 ,C 
AIK NXK ) =FEAL(01 (2) ) 
A22(NXK ) = FEAL(01 ( 1) ) 

K = K + 1 

11 CONTINUE 
1C CO 12 I=K,NCIF 
N>K = NX1< + 1 
ACCINXK ) = 1.C 
AIKNXK )=C ,C 
A22(NXK ) = C,C 
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11 CC.N7IME 
GC TC 2C 
3 IX = l 

CG 13 1 = 1, ^XP 
AC0( I ) =1 ,C 
;ill( I) = XF (I) 

^22( I ) =C .C 
1* CGNTINL5 
NXX = i\XF4l 

CC 14 i=N>>,yx 
I>l=IX4l 

C(1)=C^FLX<VS(IX^,Y!(IX)) 
C(2)=C^FL> (Ys(lXi> ,YI(IX1) ) 
CALL MAKFCL (2,0, Cl) 

ACC(NXY ) = 1 ,C 

AIK NX> ) =?6AL(D1 (2) ) 

A22(NX> )=FEAL(Di( 1) ) 

IX = IXH1 
N>X = NX>+ 1 
14 CCN7I.NLE 
2 C FE 7U F N 
END 



SL3RCLTINE PCL£XP(XX , 1 0 1 1"* X X , M , YY , I C I ^ Y Y ) 

SL3PGL7INE TC FIND 7HE FCwEFS CF GIVEN FCLYNOV I^L 

X> : FCLY.NCMAL 

ICIVXX : CI^c^SICN CF XX 

y: PQV^ER TC EE RAISED 

YY; RESULTANT FCLYNGMIAL 

ICINYY: CIyE^SIC^ OF RESULTANT PCLYNGMIAL 

CINENSICN XX(25) ,YY(2E),ZZ(25) 

IF(M.EC.l) GC TO 4 

ICIMYY=1 

YY(1 ) = 1.C 

IF(y.EC.C) GC TC 3 
CC 1 I=1,Y 

CALL FNFY (ZZ,IDiyZZ,YY,ICiyYY,XX,IDI‘‘‘XX) 

CG 2 j=i,iciyzz 

YY( J l = 2Z( J ) 

2 CCNTIME 

ICiyYY=ICI>ZZ 
1 CGNTINLE 
GC TC 2 

4 CQ 5 I = 1,ICKXX 
YY ( I l = XX ( I ) 

E CCNTINLE 

ICIMYY = IC lYXX 
I FETUFN 
END 



SL3RCLTINB PPLGT (RC , RC 1 , W B EGI N , « CELT A , N cC ! NT , PX , P Y , 
=^NCi,yN n 

DIVENSICN FC (100) ,RC1 (ICO) 

CINENSICN FX(IJO) ,?Y(10C) 

CCyPLEX XMy,CNUM 
'*» = ^8EG IN 

CC 666 I=1,NPCINT 
> = '*^ + VC£LTA 

INITIALIZATION CF REAL ANC CC^PLEX PARTS JF NU'^EPATOR 
AND DENCMNATCF 

XNUMA = C ,C 
XNLye = C. C 
XNLyC = C .C 
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XNLMD=C.C 

XCNUNA=C .C 
X0NUNB=C.C 
XDNUMC=C .C 
XCNUVC = C .C 

IMTI4LIZ4TICN OF CQLMEPS 

^NLMC=1 
NNUM 1=2 
NNUM2=2 
MNLM3=4 

/'CNUNC = 1 
NDNUM1=2 
FCNUN2=2 
N0NUN3=4 

XG=MM.>C-1 
X1=M^L^ 1-1 
X2 = MNLf'2-l 
X3 = f^M.>3-l 

XC0 = ^CMFC-1 
XCl = ^C^LM-l 
XD2 = ^D^^.f'2-l 
XC3 = f'C^U^3- 1 

00 6£5 J = 1 ,100 

IF(MMF3 .LE .^<N1 ) XNUNC=XNUNC+RC ( MNUM3 ) =>W**X3 
IF(MNUf< 2 .LE .MNl) XNL NO = XNL VC + RC ( MNUM2 ) =^W«*X2 
IF(MhUM .LE ..‘^Nl ) XNUf'B=XNUNB+RC(MNUM 1 
IF(MNLNC.LE.NNl) XNUN A = X^U^'A + RC (MNUMC )=5'W’*«XC 

IF (MCNL^B .LE .MOl ) XDNUMC=XCNUMO + RC U MDNLM3 ) ’^V»-=5‘X03 
IF(MCMV2.L£.MOi) XC^L^C = X C^UMC+RC1 ( VCNUM2 ) =»W=S‘>i=X 02 
IF(yCNLMl .LE.MOl ) XOMjNB=XCNUMB + RCl ( )=»V»*=!'XD1 
IF( yCMVC .LE .MOi ) XC^L■M A = X CNUMA + RC 1 ( MCNiy 0 ) =» U*-X00 

INCRSyEM COUNTERS 

NMHC = NNLNC + ^ 

NNUMl = NNUM+4 
NNUy2=NNLN2+4 
MNUM3=NNLN3+4 

NDNUNC=NCNLyO+4 
yONUMl = yCNLM + 4 
yCNUy2=NCNLy2+4 
N0NUN3 = >-CNLN3+4 

X0 = XC .0 
X1=X1+^.C 
X2=X2+4 .C 
X3=X344.C 



XC0=XCC + ^ .C 
XCl = XCl-»4 .C 
>C2=XC2+4.C 
XC3 = XC3 + ^ .C 

IF( (NDNLNC.GT.MOl) .ANC. (MNLyO.GT.MNl ) ) GO TO 668 
66: CONTINUE 

CAuCUiLATE NAGMTUDE OF 7FE TRANSFER FUNCTION 

66E XNUy 1=>NUNA-XNUMC 
XNLM2=>NUNE-XNUM0 
XNUM = CNFLX ( XNUM 1 , X NU N 2 ) 

CNUN l=>ONLNA-XONUyC 
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C^UM2=>C^L^E-X0NUM0 
CNU^' = C^PL> (CNUNl ,CNUN2) 
>MAGh=CAE« (XNUM/ONLM ) 
FX(I \ = V 

PVd )=2C.C>>AL0Gl0(xyAGH) 
66^ CONTIME 
FETUFN 
END 



SCaRGLlINE XFQRM(X,Y,MYPEtXT,VT,M,N,Fl,F2 ) 
SL8RGUTINE FCF FREQUENCY TRANSFGRyATIONS 



N7YFE=0 
NTYPE=i 
N1YPE=2 
MYP6 = 3 
N1YPE=4 
X 
Y 
M 
N 

XT, YT 



NG 7RANSF0RMATICNS REQUIRED 
LCV»-PASS TO LCVi-PASS TRANSFORMAT ICN 
LCl^-PASS TG HIGH-PASS TRANS FO RN AT ION 
LCV'-PASS TO BANC-PASS TRANSFORMATION 
LCy«-FASS TC EANC-STCP TRANSFORMAT IGN 
COEFFICIENTS OF NUMERATOR PCLYNCMIAL 
CCEFFICIENTS OF CENCMINATGR POLYNOMIAL 
ORDER OF NLMERATCR PCLYNCMAL 
ORDER GF CENOMINATOR POLYNOMIAL 

TRANSFORMED NUMERATOR AND DENOMINATOR PCLYNOMIALS 



CIMENSICN X(25),Y(25 »tXT(25),YT(25) 



NM=N-M 

N1=N+1 

M1=M+1 

IF(NTYFE.NE.l) GG TC 11 
CO 1 1 = 1, Ml 
XM=N l-I 

XT(I ) = > (I ) =» (F1**XM) 

1 CCNTINLE 

CC 2 1 = 1, N1 
>N=N1-I 

YT( I ( I d (F1^==*XN ) 

2 CCNTINLE 
GO TO 2C 

11 IF(NTYFE .NE .2 ) GO TO 20 



XF=M 

YF=N 

FACT = ( Fl=»=*\F)/(Fl«=i‘XF) 

K=M1+1 

CG 3 1 = 1, Ml 

XM=M1-1 

K = K — 1 

XT (I > = (X (K »4(F1*=«=XM ) )*FACT 
3 CCNTINLE 
K=N1+1 
CC A 1 = 1, N1 
XN=N l-I 
K = K 1 

YT(I ) = 'r (K )4(F1**XN ) 

A CCNTINLE 
2C RETURN 
END 



SUBROUTINE ROOT ( N , RR ,R I I 

SL3R0LTINE TO FIND THE PCCTS CF THE EQUATION 
X«^N+1=0 

RF* lARRAY CC NT A INI NG ^ RE AL PARTS OF CALCLLATED PCCTS 
RI :ARRAY CONTAINING IMAGINARY PARTS OF CALCULATED ROOTS 
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CI^E^SK^ FR(25» ,RI(25) 

PI=3. l^lf ?2 

XN=N 

T5TA=FI/XN 
OC 1 1=1, N 
XI=I-l 

ARG=( (2.C=»XI ) + l)«TETA 
RR( I ) = «U (AFG) 

RI (I > = CCS (ARG ) 

1 CCMIME 
RETURN 
END 
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APPENDIX B 



COMPUTER PROGRAM LISTING 



COMMCN H(ICC) ,HM(lC0,100)tX(100) ,V ( 1 00 ) ,N, P ,NPULS E 

C SLoRCUTINE TC FIND THE TIME DOMAIN RESPONSE 
C A DIGITAL FILTER 
C 

C N : FILTEF'S CRCER + 1 

C A. : NUMERATOR OOEFFICIEMS OF THE DIGITAL FILTER 
C e ; DENOMINATOR COEFFICIENTS ( WHERE FIRST COEFFICIENT 
0 WIuL BE NORMALIZED AND WILL NOT EE ENTERED TC THE 
C PROGRAM 

0 

F = l. 

N = 24 

NFULSE=2C 
WRITE , 1C) 

1C FQRMAT('l') 

CALL IMFLLS 
CALL HMTRX 
CALL INFLT 
CALL CCNVCL 
WRITE ( t , 1C) 

CALL FLCT 
WRITE(^ ,1C) 

STOP 

END 



SLERCLTINE HMTRX 

SLBROUTINE TC FIND SYSTEMS TRANSFER MATRIX 

COMMCN H (100) ,HM( 100 ,100) ,X(100) ,V(100) ,N,P,NPULSE 



NN=N+1 
NM=NFLLSE+1 
DC 1 1 = 1, NN 
DC 2 J = 1 ,NN 
H,M(I ,J)=C.C 
2 CONTINUE 
1 CONTINUE 
CO 2 1 = 1, NN 
CO 4 J = 1 ,NN 
IF(J.G1.I ) GO TO ^ 
L=I-J+1 



HM(i ,L ) = H(J) 

4 CONTINUE 
3 CONTINUE 

WRIT£(C,2C ) 

2C F0RMAT(2X,' HMATRIX •///) 

3C FORMAT (2X ,26F5. 3// ) 

CC 5 1 = 1, NN 

WRITE(6,3C) (HM( I , J) ,J = 1 ,NN ) 

5 CONTINUE 
RETURN 
END 



SLBRCLIINE INPUT 

SLBROUTINE TC FORM THE INFLT PULSE SEQUENCE 

COMMON H(1C0),HM(100,100),X(10C),V(1C0) ,N, P ,NPULSE 
NN=N+1 
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NM=NFULSE+1 
CC 1 I=1,NFLLSE 
>( I ) = F 

1 CCNTINLE 
CC 2 I =hN 
X(I )=0.C 

2 CCNTIM.E 
VRITE ( 6 t 1C ) 

1C FCRMAT(2X,' INPUT VECTOR ’//) 
V»RITE(£,2C) (X(I ) ,I = 1,NN) 

2C FORMAT (2> ,26F5. 3// ) 

RETURN 

END 



SLERCLTINE CCNVOL 
C 

C SL3ROUTINE TC PERFORM CONVOLUTION SUMMATION 
C 

COMMON F (100) ,HM( 100 ,100) , X(130) ,V (100) ,N,P,NPULSE 

NN=N*1 

CO 1 1=1, 

V(I )=0.C 
CG 2 J=l,^^ 

V(I )=V(I ) + l-N(I ,J)’«‘X(J) 

2 CCNTINLE 
1 CCNTINLE 
VnRITE(6 ,1C) 

10 F0RMAT(2X,' RESULT VECTOR *//) 

)sRITE(A,2C) (V(l ) ,I = 1,NN) 

2C FORMAT (2> ,2£F5. 3//) 

RETURN 

ENC 



slercltine plot 

C 

C SLBRCUTINE TC FLQT THE OUTPUT 
C 

COMMON h ( 100) ,HM( 100 ,100) , X( 100) ,V (100) ,N,P ,NPULSE 

CIM5NSICN T (250) ,Y(25C) ,Z( 250 

NN=N+1 

NCELTA=1C 

NFCINT=(NFLLSE+3)=»NCELTA 

CO 1 i=l,NFCINT 

XI = I 

T(I) =XI 

Y( I )=C ,C 

Z(I )=C ,C 

1 CCNTINLE 
K - C 

CO 2 1 = 1 ,NFCINT,NDELTA 
I F(K.GT .NN ) GC TQ 5 
K = K+1 

CG 3 J=l,4 
L=I+J-1 
Y (L ) =V (K ) 

3 CCNTINLE 

2 CCNTINLE 

5 CALL FLCTF (T,Y,NPCINT, 1 ) 

C 

C PLOT INFLT FLLSES 
C 

K = 0 

NM=NFLLSE-1 

CO 10 1 = 1 ,NPCINT ,NDELTA 
I F(K .GT .NM ) GC TG 7 
K = K + 1 

CC 8 J=l,^ 

L=I+J-1 
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Z(L>=X (K ) 

6 CGMIME 
1C CONTINLE 

7 CALL PLCTF(T,Z,NPCINT,2) 
RETURN 

END 



SLBRCLTINE INPULS 



SISROUTINE 7C FIND INFULSE RESPONSE OF THE 
( SOLVING DIFFERENCE ECLATIONS ) 



CgWMCN _K( ICC) ,HM(1C0 



CINENS ICN 

CINENSICN 

CIMENSICN 

READ (5 ,2C 

NN=N+1 

NNN=NN42C 



.LLi in 

A ( 100 ) . _ . . 
CLNMY( 100) 
C ( 100) 

) N 



. ilOO) 
B ( KC ) 



X(IOO) iV(lCO) 



FILTER. 

N, PiNPULSE 



INITIALIZATION OF VECTORS 

CO 35 1 = 1, NNN 

C(I)=C, 

h(I )=C. 

A(I )=C. 

E(I)=C. 

35 CONTINLE 
C( 20) =1 . 



I 



SHIFT ORIGIN TO 20 
NQR0ER = f- + 2C 

R£AC(5,2C)(A(II,I=20,NCRCER) 

READ (5,2C)(E(I),I=2C,NCRCEF) 

3C FCRMAT(I2) 

2C FCRNAT (8F1C.5) 

VRITE(A ,1) N 

1 FGRMAT(2)(,' ORDER OF THE FILTER IS :M2///) 
WRITE( E ,2 ) 

2 FORMAT (2> 1 'NLNERATOR COEFFICIENTS '///) 
WRIT£(e,3)(A(n,I = 20,NCRCER) 

WRITE(E iA) 

4 FORMAT (2X 1 'CENOMINATCR COEFFICIENTS '///) 
V»RITE(E,3)(E(I), 1=20 iNCRCEP ) 

CO 4C II=2C,NNN 
FACT1 = C . 

NK=I I+l 

CO 5C IL=2C,N0RD£R 
NK=NK- 1 

FACT1=FACT1+A(IL)*C(NK) 

5C CONTINLE 
FACT2 = C . 

NK = I I 

CO 6C IL=21iN0RD£R 
NK=NK-1 

FACT2 = FACT2+B( ID =Hh(NK) 

6C CONTINLE 

h( II ) =FACT 1-FACT2 
4C CONTINLE 

H I T E ( E 1 2 1 ) 

21 FORMAT (2X 1 'WEIGHTING SECUENCE * /// ) 

WRIT E( 6 1 1C KH( I ) , I = 2C,NNN ) 

1C FORMAT (2> 1 12F10. 5// ) 

2 F0R,mAT(2>i12F10.3) 

K K K = 1 

CO 55 KK=2C,N0RDEP 
A (KKK ) =A ( KK ) 
e (KKK ) =E (KK ) 

KKK = KKK + 1 
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55 CONTINLE 
K K — 1 

DO~65 I=2C,NNN 
CUMI^Y { KK )=l- ( I » 

KK=KK+1 
6! CGNTINLE 

CO 75 1 = 1, NN 
hdJ =CLNNY (I ) 

75 CONTINLE 
>CUM = e (2 ) 

B(3) =B (1 J 
e(l)=XCLN 
READ (5 , ICC J F1,F2 
P£A0(5,1CC) T 
IOC FORMAT (2F1C. 5) 

NF=1CC 
>NP=1CC . 

FCELTA= ( F2-F1 )/XNP 

CALL 2CMN(A ,E ,M , M , F 1 , F2 , FDcLT A ,N P , T ) 

RETURN 

ENC 
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